Mon. Not. R. Astron. Soc. 000.11111(31(20061 Printed 5 February 2008 (MN WF&t style file v2.2) 



A Subaru/ Suprime- Cam wide-field survey of globular cluster 
populations around M87 - I: Observation, data analysis, 
and luminosity function 



v,o ;Naoyuki Tamura 1 *"]", Ray M. Sharpies 1 , Nobuo Arimoto 2 , Masato Onodera 2 ' 3 , 



Kouji Ohta 4 , & Yoshihiko Yamada 2 



O 

o 

^v^j ■ 1 Department of Physics, University of Durham, South Road, Durham, DH1 3LE, United Kingdom 
2 National Astronomical Observatory of Japan, Mitaka, Tokyo 181-8588, Japan 
Or 3 Department of Astronomy, School of Science, University of Tokyo, Tokyo 113-0033, Japan 
Department of Astronomy, Faculty of Science, Kyoto University, Kyoto 606-8502, Japan 

C/3 , 



^ ! ABSTRACT 

\o '. 

f- ^ In this paper and a companion paper, we report on a wide-field imaging survey 

■ of the globular cluster (GC) populations around M87 carried out with Suprime-Cam 
' on the 8.2m Subaru telescope. Here we describe the observations, data reduction, 
, and data analysis and present luminosity functions of GC populations around M87 

■ and NGC 4552, another luminous Virgo elliptical in our survey field. The imaging 
"^h \ data were taken in the B 7 V, and / bands with a sky coverage of 2° x 0°5 extending 

from the M87 centre out to ~ 0.5 Mpc. GC candidates were selected by applying a 
I ' colour criterion on the B — V and V — I diagram to unresolved objects, which greatly 

reduces contamination. The data from control fields taken with Subaru/Suprime-Cam 
■£j i were also analyzed for subtraction of contamination in the GC sample. These control 

field data are compatible with those in the M87 field in terms of the filter set (BVI), 
. . , limiting magnitudes, and image quality, which minimizes the possibility of introducing 

any systematic errors into the subtractive correction. We investigate GC luminosity 
J^j ' functions (GCLFs) at distances ^ 10' (< 45 kpc) from the host galaxy centre in 

detail. By fitting Gaussians to the GCLFs, the V-band turnover magnitude (Vto) is 

■ estimated to be 23.62 ± 0.06 mag and 23.56 ± 0.20 mag for the GC population in 

M87 and NGC 4552, respectively. The GCLF is found to be a function of GC colour; 
Vto of the red GC subpopulation (V — I > 1.1) is fainter than that of the blue GC 
subpopulation (V — I ^ 1.1) in both M87 and NGC 4552, as expected if the colour 
differences are primarily due to a metallicity effect and the mass functions of the two 
subpopulations are similar. The radial dependence of the GCLF is also investigated 
for the GC population in M87. The GCLF of each subpopulation at V ^ R ^ 5' is 
compared with that at 5' ^ R ^ 10' but no significant trend with distance is found in 
the shape of the GCLF. We also estimate GC specific frequencies (Sn) for M87 and 
NGC 4552. The S N of the M87 GC population is estimated to be 12.5 ± 0.8 within 
25'. The S N value of the NGC 4552 GC population is estimated to be 5.0 ± 0.6 within 
10'. 
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Globular clusters (GCs) are homogeneous stellar systems 
containing stars with a single age and metallicity, which 
are in principle simpler to interpret than photometric and 
spectroscopic observations of the integrated stellar light of a 
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galaxy. GCs are therefore considered to be powerful probes 
with which to understand the star formation and chemical 
enrichment history of their host galaxy. One of the basic 
findings from observations about GC populations in lumi- 
nous galaxies is that while thousands of GCs are associ- 
ated with luminous elliptical galaxies, a significantly smaller 
number of GCs exist around spiral galaxies with similar lu- 
minosities (e.g., Harris 1991; Barmby 2003). This indicates 
that the specific frequency of GCs (Sn), which is consid- 
ered to be related to the relative efficiency of GC formation 
and/or survival compared to galactic halo/bulge stars, de- 
pends on galaxy morphology. In fact, Sn has also been sug- 
gested to be correlated with local galaxy density, with galax- 
ies in denser environments having larger Sn values (West 
1993). The fact that this trend appears to exist even when 
the sample of galaxies is restricted to ellipticals suggests 
that GC formation efficiency is more physically linked with 
galaxy environment. One possibility to explain this observa- 
tion is biased GC formation in galaxies inhabiting denser en- 
vironments (West 1993; Blakeslee 1999; McLaughlin 1999). 
Alternatively, a substantial number of GCs in a luminous 
galaxy may have an external origin; GCs could be captured 
from other galaxies through galaxy interactions and accrete 
onto luminous galaxies in clusters, which would then en- 
hance Sn values as observed. 

Recent studies of GCs in the central regions of lumi- 
nous ellipticals conducted with the Hubble Space Telescope 
(HST) have revealed that many luminous ellipticals have 
bimodal or multimodal colour distributions of GCs (Geb- 
hardt & Kissler-Patig 1999; Larsen et al. 2001; Brodie et al. 
2005). It is found that the mean colours of both red (metal- 
rich) and blue (metal-poor) GC subpopulations are corre- 
lated with the host galaxy luminosities and colours (Larsen 
et al. 2001; Strader, Brodie & Forbes 2004; Strader et al. 
2005; Peng et al. 2006), which will be important constraints 
on the proposed scenarios for the formation and evolution of 
GC population such as multiphase collapse scenario (Forbes, 
Brodie & Grillmair 1997), merger scenario (Ashman & Zepf 
1992), hiearchical merging scenario (Beasley et al. 2002), 
and accretion scenario (Cote, Marzke & West 1998). The 
small field of view of the HST, however, does not allow one 
to collect GC populations in the outer halo of a galaxy and 
investigate their spatial structures, which will also be key 
pieces of the puzzle (e.g., Moore et al. 2005). Much wider- 
field studies (e.g., out to ~ 100 kpc from the galaxy cen- 
tre) of GC populations therefore need to be performed us- 
ing ground-based data (Rhode & Zepf 2001, 2004; Dirsch et 
al. 2003; Bassino et al. 2006). Studying GC populations at 
large distances from the host galaxy is of great importance 
because the outer halo of the host galaxy, and even the in- 
tergalactic space, are presumed to be large reservoirs of blue 
and metal-poor GCs in the accretion scenario. 

In this paper and a companion paper (Tamura et al. 
2006; Paper II hereafter), we report a wide-field imag- 
ing survey of the GC populations around M87 conducted 
with Suprime-Cam on the Subaru telescope. Several mod- 
erately wide-field studies of GCs surrounding M87 have 
already been carried out using photometry (Strom et al. 
1981; McLaughlin, Harris, & Harris 1994; Harris, Harris, 
& McLaughrin 1998; Hanes et al. 2001) and spectroscopy 
(Cohen & Rizhov 1997; Cohen, Blakeslee & Rizhov 1998; 
Kissler-Patig & Gebhardt 1998; Cote et al. 2001), but all of 



these studies explored only the regions < 10' (< 45 kpc or 5 
r e ) from M87. In contrast, the area of our survey is approx- 
imately 2° x 0°5 (560 kpc x 140 kpc) extending from the 
centre of M87 out to ~ 0.5 Mpc, which is the widest survey 
yet undertaken of the GC populations in luminous galaxies. 
In this paper, we focus on describing the observations, data 
reduction (§|5J| and data analyses such as selection and pho- 
tometry of GC candidates, incompleteness correction, and 
subtraction of foreground and background contamination in 
the GC candidates (§EJ. We derive GC luminosity functions 
around M87 and NGC 4552, another Virgo luminous ellip- 
tical galaxy in our survey field, and estimate the global GC 
specific frequencies of these luminous ellipticals in § 2] We 
investigate colour distributions and spatial distributions of 
GC candidates in Paper II. We adopt distances of 16.1 Mpc 
(distance modulus of 31.03) to M87, and 15.4 Mpc (distance 
modulus of 30.93) to NGC 4552, based on measurements us- 
ing the surface brightness fluctuation method (Tonry et al. 
2001). An angular scale of 1' corresponds to 4.7 kpc and 4.5 
kpc at the distance of M87 and NGC 4552, respectively. 

2 OBSERVATIONS, DATA REDUCTIONS, 
AND CALIBRATIONS 

2.1 The M87 Fields 

2.1.1 Observations and data reductions 

Imaging observations were performed on 17 and 18 March 
2004, with Suprime-Cam (SCam) (Miyazaki et al. 2002) on 
the Subaru telescope. SCam is a mosaic CCD camera with 
10 2Kx4K CCD chips and the field of view is approximately 
34' x 27' on the sky. The pixel scale is 0''2 pixel" 1 . In this 
observing program, a field of approximately 1 square degree 
(136' x 27') extending from M87 towards the east was cov- 
ered by 4 telescope pointings through B-, V-, and /-band 
filters. The field IDs and locations are shown in Fig. and 
the observation log is presented in Table Each field was 
observed with the telescope dithered by ~ 5". Since this 
dithering scale is smaller than the gap between the CCDs, 
the 10 CCD frames are not mosaiced into one continuous 
frame but are reduced and analyzed individually. A typical 
exposure time of one frame is 360 sec, 270 sec, and 240 sec 
in B, V, and I band, respectively; several frames were co- 
added to give the total exposure times listed in Tabled On 
the first night, the sky condition was non-photometric and 
the transparency was highly variable. On the second night, 
it was much better but was still hazy with a little varia- 
tion. We therefore scale the data taken on the first night 
by shifting the magnitude zeropoints to match with those 
of the data on the second night, calibrate the reduced data 
based on the standard stars taken on the second night, and 
check the calibration using GC photometry in the literature 
(see next section for details). Typical seeing sizes during the 
observations were l"5 in B band and ~ l"0 in V and / 
bands. 

Data reduction was performed with IRAF 1 in a stan- 

1 IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which is operated by the Association of Universities for 
Research in Astronomy, Inc. under cooperative agreement with 
the National Science Foundation. 
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Field ID Field centre Integration time (mii m ) Seeing size 



[sec (mag)] [arcscc] 



(1) 




(2) 






(3) 






(4) 






a(J2000) 


<5(J2000) 


b(J2000) 


B 


V 


I 


B 


V 


I 


Field 1 


12 h 31 m 18 s 4 


12°29'13" 


74° 6 


3680 (25.6) 


1350 (25.1) 


3480 (24.6) 


1.8 


1.0 


1.0 


Field 2 


12 h 33 m 42 s 8 


12°27'37" 


74° 8 


2640 (25.9) 


1350 (25.2) 


3690 (24.7) 


1.2 


1.0 


1.0 


Field 3 


12 h 36 m 08f8 


12°24'35" 


74° 9 


1800 (25.7) 


1350 (25.2) 


2640 (24.5) 


1.5 


1.0 


1.1 


Field 4 


12 h 38 m 33 s 2 


12°24'35" 


75°0 


2280 (26.1) 


1350 (25.2) 


4380 (24.8) 


1.1 


1.0 


1.0 


HDF-N 
Lockman Hole 


12 h 36 m 46 s 7 
10 ,l 35 m 55f2 


62°11'50" 
57° 42' 18" 


54° 8 
51°3 


6000 (26.6) 
6000 (26.9) 


4800 (25.5) 
4800 (26.0) 


4200 (25.2) 
3600 (24.8) 


0.8 
1.0 


1.1 
1.1 


0.9 
1.3 



Table 1. Observation log. Data for HDF-N and Lockman Hole are retrieved from SMOKA. Col. (3): The integration 
times were calculated including data taken under non-photometric conditions on the first night (hence the interpre- 
tation is not straightforward). Number in the parentheses indicates 50 % completeness to point sources estimated 
with artificial star test. Col. (4): The seeing sizes are estimated in the stacked images. 



.- 



i> m S 

■ 

I m ' — * '. 



- i.' ■ 



d 4 Held 3 

i ■ <*n~- — i — — . 



J 



NCC 4552 



Held. 1 - 
1 



: 



Figure 1. Locations and IDs of the observed fields. Each box indicates one SCam field of view. The background is a DSS image. North 
is up and East is left. 



dard manner; bias subtraction, fiat-fielding, masking bad 
columns and saturated pixels, sky subtraction, registration, 
and average stacking with a 3 a clipping algorithm. Sky 
subtraction was performed by employing the following two 
steps. Firstly, an image was divided into a mesh of 128 x 
128 pixels (~ 26" x 26") and a median sky value was esti- 
mated in each window after bright objects were masked. A 
sky value in each pixel was then estimated by an interpola- 
tion of the median sky values for the adjacent windows. A 
background image of an object frame was created with this 
process which was then subtracted. For CCD frames where 



bright galaxies or their envelopes are quite extended (e.g., 
near M87 and NCC 4552), this method cannot be applied. 
Instead, an average background was estimated as a single 
value using a "blank" CCD frame within the same expo- 
sure. It was corrected for the sensitivity difference between 
the two CCD frames using the flat-field frames. In stacking 
CCD frames, aperture photometry of 20 — 30 bright stellar 
objects selected using SExtractor (Bertin & Arnout 1996) 
based on the CLASS_STAR index was performed in each frame 
and the zeropoint of the frame was shifted so as to match 
that in a frame taken on the second night. PSF matching 




Figure 2. Reduced data for Field 1 (V band). North is up and 
East is left. 



Figure 4. Same as Fig. but for Field 3. 
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Figure 3. Same as Fig. [5] but for Field 2. 

was not performed to avoid degradation of image quality. 
The stacked images are presented in Fig. [5] 

2.1.2 Photometric and astrometric calibration 

Photometric calibration of the M87 fields was performed 
using standard stars from Landolt (1992) which were ob- 
served at the beginning and end of the second night. Several 
standard stars were imaged on each CCD chip so that the 
calibration could be carried out individually. Photometry of 
standard stars was performed within a 12" diameter aper- 
ture. After excluding saturated stars and those in crowded 
regions, magnitude zeropoints and the trends with sec z and 
colours were calculated. The estimated accuracy in the fit- 
ting procedure is 0.03 — 0.05 mag. 

Since our data include frames which were taken under 
non-photometric sky conditions, we check the photometric 
calibration by comparing our photometry of GCs around 
M87 with that in the literature. We make use of GC photom- 
etry in the V and I bands by Harris, Harris, & McLaughlin 
(1998, HHM hereafter) and that in the Washington system 
(C and Ti) by Hanes et al. (2001, H01 hereafter). The GC 
photometry by H01 is converted to the standard BVI sys- 
tem by using the formulae obtained by Geisler (1996). PSF- 



Figure 5. Same as Fig. [5] but for Field 4. 

fitting photometry was carried out for GCs using the PSF 
and ALLSTAR tasks in the DAOPHOT package of IRAF. 
A PSF was determined using ~ 20 moderately bright (un- 
saturated) stellar objects in each CCD frame, which were 
selected using SExtractor based on the CLASS_STAR index. 
The PSF obtained was also used for GC photometry, since 
the GCs at the distance of the Virgo cluster are unresolved at 
the seeing sizes of our images. Galactic extinction was then 
corrected using reddening maps from Schlegel, Finkbeiner, 
& Davis (1998). We are primarily concerned about zeropoint 
offsets for l/-band magnitude and B — V and V — I colours; 
the V band image is used as a selection band to make a cat- 
alog of GCs and B — V and V — I colours are used to isolate 
GC candidates from other unresolved objects (see § 13.31 for 
details). While C — Ti can be converted into B — V with 
a small error, it is not converted to V — I with a good ac- 
curacy (Geisler 1996). We therefore decided to estimate the 
zeropoint offset in B — V from the comparison with H01 
and those in V and V — I from the comparison with HHM. 
Our GC photometry is compared with that from H01 and 
HHM in Fig.|S| Dotted lines in these plots indicate the ap- 
proximate edges of the colour range which is expected to be 
occupied by GCs. These comparisons suggest that there are 
some zeropoint offsets between our photometry and that in 
the literature. The zeropoint offsets in V, B — V and V — I 
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Figure 6. Comparison of our photometry (VsCam, (B — V)sCami and (V — -OsCam) with that by HOI ((B — V)hoi) and HHM (Vhhm 
and (V — /)hhm)- Dotted lines in the plots for the colours approximately show reddest and bluest boundaries expected for GCs. 



are 0.04 mag, 0.05 mag and 0.10 mag, respectively. Our mag- 
nitude and colours are corrected for these zeropoint offsets 
in the following analyses. 

We performed astrometric calibration against the 
2MASS catalog using the CCMAP task in IRAF. A plate 
solution (second order polynomial with full cross terms) was 
computed for each CCD chip using stars with J ^ 17 mag. 
The fitting accuracy is typically ~ 0"2. 



2.2 The Control Fields 

In addition to the M87 field data, we also analyze BVI im- 
ages of the (blank) fields, HDF-N and Lockman Hole (LH), 
retrieved from the Subaru-Mitaka-Okayama-Kiso Archive 
(SMOKA) system (Baba et al. 2002). These data were taken 
with SCam on the Subaru telescope during several observ- 
ing runs in 2001: 23 and 24 Feb for BVI on HDF-N and 
B on LH, and 22 and 23 Apr for V and / on LH (Capak 
et al. 2004). This information is also summarized in Table 
These control field data are used to estimate contami- 
nation in the sample of GC candidates and to statistically 
subtract it, which is essential for investigating properties of 
GC populations such as the luminosity function and colour 
distribution. We emphasize that these control fields are also 
at high galactic latitudes (Table 0, and that the data cover 
reasonably wide sky areas (one SCam field of view: ~ 900 
arcmin 2 ) and are comparable with our data in the M87 field 
in terms of the filter set, limiting magnitudes, and image 
quality. This minimizes the possibility of introducing any 
systematic errors into the subtractive corrections for fore- 
ground and background contamination in the GC sample. 

The data reduction was carried out with SDFRED 
(Yagi et al. 2002; Ouchi et al. 2004), which is a reduction 
pipeline optimized for SCam data of blank fields. The ba- 
sic reduction procedure is the same as that applied to the 
M87 field data. The large scale of telescope dithering (~ 1') 



for these control field data enables the CCD frames to be 
stacked into one continuous image with sensitivity differ- 
ences between CCDs corrected by using stellar objects in the 
overlap regions. To avoid any complications due to a possible 
drop of limiting magnitude near the field edge, we use only 
the central 27' x 27' region in the following analyses. The 
magnitude zeropoints of the HDF-N data were calculated 
using the photometry catalog by Capak et al. (2004) . These 
authors did not observe any standard stars during the ob- 
servations and determined magnitude zeropoints by exploit- 
ing the accurate photometry of objects in the region where 
deep HST/WFPC2 data are available. The best-fit SEDs to 
the multi-band photometry of the objects (Fernandez-Soto, 
Lanzetta & Yahil 1999) were used to account for the slight 
differences in the filter responses between Subaru/SCam and 
HST/WFPC2. Since the B band data of the LH field were 
taken in the same runs as for the HDF-N, the same zeropoint 
is adopted. The V and / band data of the LH field were 
taken on different observing runs and standard stars were 
observed at elevations similar to those of the LH field dur- 
ing the night. Magnitude zeropoints were derived from these 
data. Galactic extinction was then corrected using Schlegel 
et al. (1998). 



3 DATA ANALYSES 

3.1 Halo Light Subtraction 

As shown in Figs. H and |1] M87 and NGC 4552 extend 
across several CCD frames and their halos have to be sub- 
tracted to reveal the GC populations. We removed the halos 
by conducting an iterative median smoothing and subtrac- 
tion (e.g., McLaughlin, Harris, & Hanes 1994). First, we 
subtract unresolved objects from an image. This process is 
not mandatory, but it helps better model the extended halo 
light distributions or their residuals in subsequent iterations. 
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Figure 7. Overall detection completenesses in BVI bands are overplotted. Top and middle panels: the results in the M87 fields are 
shown. Bottom panels: those in the control fields are indicated. 



Object detection was performed with SExtractor and unre- 
solved objects were picked out based on the CLASS_STAR in- 
dex. A PSF was determined using moderately bright stars, 
which was then fitted to unresolved objects and the fitted 
profiles were subtracted from the original image using the 
ALLSTAR task. The resulting image was median-smoothed 
to create an image with a halo light distribution which was 
subtracted from the original image. This procedure was re- 
peated 4 times with the mesh size successively reduced (f 28, 
64, 32, and 32 pixels). 

3.2 Artificial Star Tests 

We performed artificial star tests for investigating detection 
completeness to point sources on the B-, V-, and /-band 
images, after bright galaxies were subtracted if necessary. 
Using the fRAF STARLIST and MKOBJECT tasks, artifi- 
cial stars with the same PSF as that determined using real 
unresolved objects were distributed on the original image. 
For a series of tests, 500 artificial stars within a certain range 
of magnitude [m, m + 0.5 mag] were generated while succes- 
sively changing the magnitude range as the test progressed. 
SExtractor was then used for object detection and the frac- 
tion of artificial stars detected (i.e., detection completeness) 
was calculated as a function of magnitude. Note that faint 
stars could be rejected by DAOPHOT even if they are de- 



tected by SExtractor, but the fraction of such artificial stars 
turns out to be small f %) throughout the magnitude 
range investigated in our artificial star tests. We performed 
these artificial star tests on all the CCD frames in all the 
observed fields (Field 1 — 4) and the overall completeness 
in each observing field is plotted against B-, V-, and /-band 
magnitude in the top and middle panels of Fig. [7] This indi- 
cates the presence of a slight field-to-field variation in limit- 
ing magnitude. The detection completeness is also a function 
of galactocentric distance; for instance, limiting magnitudes 
are ~ 0.5 mag brighter near the centre of M87 where the 
noise is higher. Therefore, when investigating the luminos- 
ity function and colour distribution of GC candidates within 
an annulus at a certain distance from the host galaxy, we 
correct for incompleteness using the completeness functions 
estimated within the same annulus. Especially near the lu- 
minous ellipticals, we divide the annulus into sub-annuli to 
follow the local variation of the incompleteness in the annu- 
lus. The artificial star tests were also executed on the control 
field images and the overall completeness functions are indi- 
cated in the bottom panels of Fig. [T] The magnitudes giving 
50 % completeness on the M87 fields and the control fields 
are listed in Table Q 

Since unresolved objects are firstly selected based on the 
CLASS_STAR indices on the l/-band image when selecting GC 
candidates (see S I3.3L one also needs to consider biases asso- 
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Figure 8. CLASS_STAR indices of detected objects in the V band images of the M87 fields, HDF-N field, and LH field are plotted against 
their magnitudes in the top, middle, and bottom panels, respectively. Dashed line indicates the lower limit for selection of unresolved 
objects (0.6). 



dated with the selection of unresolved objects in addition to 
the simple detection incompleteness as mentioned above. In 
the top panels of Fig. |H] CLASS_STAR indices of the detected 
objects (left panel) and the artificial stars (right panel) in 
the M87 fields are plotted against V-band magnitude. In the 
middle and bottom panels, the results in the control fields 
are shown. These plots indicate that the CLASS_STAR index 
of a point source tends to be underestimated at fainter mag- 
nitudes and more stellar objects are expected to be excluded 
when we classify objects with CLASS_STAR indices larger than 
a certain value as unresolved. We quantify this selection ef- 
fect as functions of V magnitude and distance from the host 
galaxy based on these artificial star tests. 

3.3 Selection and Photometry of GC Candidates 

We begin selection of GC candidates with object detection 
using SExtractor. Firstly, we picked out all objects having 
at least 20 connected pixels (~ 0.8 arcsec 2 , which is approx- 
imately equal to the FWHM area of the PSF) more than 2 a 
above the local background (only objects selected in all the 
three B-, V-, and /-bands were used for analysis). Secondly, 
we selected objects with CLASS_STAR indices larger than 0.6 
as unresolved objects. The V-band image was used for this 
classification because it has the best image quality in our 



data. The subsequent results do not change if this cutoff for 
the CLASS_STAR index is set to 0.5 or 0.7. We note that our 
criterion (CLASS_STAR ^ 0.6) is more stringent than those 
adopted in previous GC studies 0.4 in Dirsch et al. 2003 
and ^ 0.35 in Forbes et al. 2004) using data with better 
image quality than our data. 

PSF-fitting photometry was performed for these unre- 
solved objects and their B, V, and I magnitudes were mea- 
sured. In this process, a residual sky background around an 
object is estimated within an annulus 10" away from the ob- 
ject with a 4" width and is subtracted. Errors in the V-band 
magnitude and B — V and V — I colours due to a PSF fitting 
error and sky subtraction error are plotted as a function of 
V-band magnitude in the top, middle, and bottom panels of 
Fig.Hfor unresolved objects in the M87 fields, HDF-N field, 
and LH field, respectively. A colour criterion was then im- 
posed on these unresolved objects to isolate GC candidates. 
We show this colour criterion as the shaded region in Fig. 
HOI which includes almost all the Galactic GCs in the catalog 
by Harris (1996) 2 but minimizes the contamination by fore- 
ground stars and background galaxies. The B — V and V — I 
colour-colour diagrams of the unresolved objects on our im- 



2 The catalog used here is the version last updated in 2003. 
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Figure 9. For the unresolved objects, errors of V-band magnitude (left), B — V colour (centre) and V — I colour (right) from the 
PSF-fitting photometry are plotted as a function of V-band magnitude. Grey line indicates the boundary below which 80 % of the 
unresolved objects are included at a given magnitude. In the top, middle, and bottom panels, the results in the M87 fields, HDF-N field, 
and LH field are shown, respectively. 
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Figure 10. The colour criterion to select GC candidates is indicated by a shaded region. In the left panel, colours of Galactic GCs in the 
catalog by Harris (2003 version) are plotted with circles, and stellar colours based on a stellar flux library by Pickles (1998) are plotted 
with asterisks. Galactic extinctions for the GC colours are corrected using the colour excess of each GC in the catalog. In the right panel, 
evolutionary tracks of the galaxy colours calculated using PEGASE v2.0 (Fioc & Rocca-Volmerange 1997) are indicated. Open squares 
are plotted at redshifts of 0, 0.5, 1.0, 1.5, and 2.0 on each track. Note that the galaxy colours at z = are plotted around the middle of 
this panel and go towards outer regions at higher redshifts. 



ages are indicated in Fig. 1111 the objects are divided into 
panels based on their distances from M87, apart from those 
in the HDF-N field, which are plotted in the bottom right 
panel for reference. Black dots indicate unresolved objects 
which pass the colour selection, while grey dots are those 
which do not satisfy the colour criterion. Some unresolved 
objects sitting outside the colour criterion are accepted as 
GC candidates by taking into account errors in the colours 
(e.g., Rhode & Zepf 2001); if colours of unresolved objects 
can satisfy the colour criterion within their errors, they are 
sampled as GC candidates. This "inclusive" colour selection 
allows us to incorporate fainter GCs which are more likely to 



be scattered out of the colour criterion due to larger errors 
in the colours. Although contaminating objects may also be 
included, they are expected to be corrected for by the control 
field data (see § 13.41 . 

3.4 Correction of Incompleteness and Foreground 
and Background Contamination 

In deriving a GC luminosity function (GCLF) and colour 
distribution, incompleteness correction is undertaken as fol- 
lows. The number of GC candidates at a certain V magni- 
tude is firstly corrected for the incompleteness in detection 
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Figure 11. B — V and V — I colour-colour diagrams for the unresolved objects. The data in the M87 fields are divided into panels based 
on the distances of the sources from M87. In the bottom right panel, the unresolved objects in the HDF-N data (see § 12.21 are plotted. 
Black dots indicate unresolved objects which pass the colour cut, and grey dots are those which do not. The shaded region in each panel 
indicates the colour criterion. Some objects sitting outside the colour criterion are accepted because we take into account errors in both 
of the colours in applying the colour cut to these objects (see text for details). 



and selection of unresolved objects on the V-band image. 
Detection incompleteness on the B and / band images is 
then corrected; GCs with a given V-band magnitude are di- 
vided into bins according to their B — V and V — I colours, 
and the numbers of objects are multiplied by a factor to 
correct for the incompleteness at the B and I magnitudes 
(B = (B - V) + V, I = V - (V - I)). These incompleteness 
corrections are also applied to GC candidates found in the 
control fields. 



Although the selection using the two colour diagram is 
expected to efficiently isolate GCs from foreground stars and 
background unresolved galaxies, there are still likely to be 
some contaminating objects, and a subtractive correction of 
this contamination is essential to investigate GC properties 
in the outer halo of the host galaxy where the GC surface 
number density is very low. We extract contaminating pop- 
ulations of unresolved objects from the control fields by us- 
ing selection criteria identical to those adopted in the M87 
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Figure 12. Upper panels: Luminosity functions of the unresolved objects which pass the colour selection in the HDF-N {left) and the LH 
{right). Thin line shows raw LF and thick line indicates incompleteness-corrected LF. Shaded region at the right hand side corresponds 
to the magnitude range where the completeness is lower than 50 %. Lower panels: Colour distributions of the unresolved objects {V $J 24 
mag) which pass the colour selection. Thin line shows raw distribution and thick line indicates incompleteness-corrected distribution. 



fields; objects with CLASS_STAR ^ 0.6 are selected as unre- 
solved objects and GC candidates are isolated by using the 
same colour criterion on the B — V vs. V — I colour-colour 
diagram. 



When subtracting contamination in a certain region 
within the M87 survey field, we take into account differ- 
ences in data quality, especially errors in the colours at a 
certain V magnitude, between the M87 field and control 
field. Since the errors in the control fields are smaller than 
those in the M87 field, a smaller number of objects would 
be scattered into the colour criterion in the control fields 
and the contamination could be underestimated. We pick 
out unresolved objects within narrow ranges (±0.1 mag) of 
V, B — V and V- 1 from both the M87 field and control field 
and calculate the differences of typical errors in the colours 
between the two samples of unresolved objects. We then ran- 
domize the measured colours of the unresolved objects in the 
control field by an amount which is determined from a Gaus- 
sian distribution whose average is zero and whose standard 
deviation is estimated from the difference of the typical er- 
rors in the colours. The observational errors in the colours 
of the unresolved objects in the control field are therefore 
replaced with the typical errors of those in the M87 field. 
This sequence is repeated in the control field for different 
V magnitudes and the colours successively changed, which 
provides a mock catalog of unresolved objects in the control 
field whose error characteristics are compatible with those in 
the M87 field. Based on the colours and errors in this mock 
catalog, the colour selection is performed to pick up contami- 



nating objects in the control field 3 . The luminosity function 
and colour distribution of these contaminating objects in- 
cluding the incompleteness corrections are then subtracted 
from those in the M87 field with the survey area normalized. 

In Fig. 1121 the V-band luminosity functions and V — I 
colour distributions of GC candidates found in the control 
fields are displayed. Note that these are obtained by con- 
sidering the entire region of the M87 survey field. The LF 
and colour distribution in the control fields are normalized 
to the survey area of a target field when subtracted. Most 
of these contaminating objects are likely to be background 
galaxies with compact morphology; candidates are dwarf el- 
lipticals and blue compact dwarfs, the latter of which need 
to be at higher redshifts than the former to meet the colour 
criterion. Foreground stars can also be scattered into the GC 
colour selection due mainly to photometric errors, but since 
the errors in the colours are ~ 0.1 mag even at V ~ 24.5 
mag, which is approximately the faintest magnitude of GCs 
studied in this work, their contribution is presumed to be 
smaller than the background galaxies. 



3 Because this sequence involves random numbers, the GC can- 
didates in a control field needs to be defined by a number of 
attempts based on the Monte Carlo technique. But in fact, the 
variance of the average population is small because the "inclu- 
sive" selection of GC candidates does not give a sharp cutoff on 
the B — V and V — I colour-colour diagram. 
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Figure 13. GC Luminosity Functions (GCLFs) obtained within annuli centered on M87 are indicated. The dotted line indicates a 
raw LF (without incompleteness correction or control field subtraction) and thin black line describes an LF after the incompleteness is 
corrected (but no field subtraction is performed yet). The red line and the thick black line are those where the control field subtraction 
is performed after the incompleteness correction: red (black) line shows the LF where the control field subtraction is performed using 
the HDF-N (LH) data, respectively. The shaded region indicates the magnitude range where the completeness is lower than 50 %. In the 
bottom right panel, the LF of unresolved and colour-selected contaminating objects in the HDF-N field is shown for reference (survey 
area is not normalized). 



4 RESULTS AND DISCUSSIONS 

4.1 Globular Cluster Luminosity Function 

In Fig. 1131 GCLFs obtained within annuli centered on M87 
are presented. The dotted lines indicate raw GCLFs with- 
out incompleteness correction or control field subtraction, 
and the thin black lines show GCLFs after the incomplete- 
ness is corrected for (but no subtractive correction for con- 



tamination is performed yet). The red (thick black) lines 
show GCLFs after the control field subtraction is also per- 
formed based on the HDF-N (LH) field, respectively. In 
Fig. 1141 the GCLFs after the incompleteness correction and 
field subtraction using the HDF-N data are divided into red 
(V — I > 1.1) and blue (V - I < 1.1) GC subpopulations. 
The black line shows the GCLF for the total GC population 
(blue + red). The boundary colour of the blue and red GCs 
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Figure 14. Black line indicates GCLF after the incompleteness is corrected and the control field is subtracted using the HDF-N data 
(same as that described with the red line in Fig, 1131 . This is divided into two LFs depending on GC colour: Red and blue lines indicate 
the LFs for red [V — J > 1.1) and blue (V — I $C 1.1) GC candidates, respectively. 



corresponds approximately to the middle of the peak colours 
in the bimodal colour distributions (see Paper II). Note that 
there is a local maximum of GC number in the 70' < r < 80' 
bin due to the GC population of NGC 4552 at a distance 
of ~ 75' from M87. GCLFs in the outermost regions tend 
to have some negative bins due to the field subtraction. In 
both of these figures, the shaded region indicates the mag- 
nitude range where the completeness is lower than 50 %. 
Note that the completeness is calculated by considering not 
only the simple detection completeness on the V band image 
but also the selection efficiency of unresolved objects based 
on the CLASS_STAR index and the completeness in B and I 



bands at the corresponding magnitudes depending on the 
GC colours. We investigate and discuss the spatial distribu- 
tions of GC populations in detail in Paper II. 

A Gaussian fitted to the GCLF obtained using only 
GCs at distances smaller than 10' (~ 45 kpc) from the 
host galaxy is plotted in Fig. 1151 (we fit a Gaussian to 
the binned data). The fainter part of the GCLF in the 
shaded region where the completeness is lower than 50 % 
is not used in this fitting process. The turnover magnitude 
(Vto) and dispersion (a) of the GCLF are then estimated 
to be Vto = 23.62 ± 0.06 mag and a = 1.50 ± 0.04 mag 
for the GC population in M87. For the NGC 4552 GCs, 
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M87 (1' < K 10') 23.62 ±0.06 1.50 ± 0.04 23.85 ±0.19 1.57 ± 0.09 23.35 ± 0.05 1.38 ±0.04 

M87 (1'<.R<5') 23.63 ±0.08 1.48 ±0.05 23.77 ±0.23 1.54 ±0.11 23.36 ±0.11 1.44 ±0.05 

M87 (5'<ii<10') 23.52 ±0.08 1.44 ±0.05 23.97 ± 0.43 1.48 ± 0.22 23.37 ±0.08 1.36 ±0.05 

M87 (iJ<l';K99) 23.67 ± 0.07 1.39 ±0.06 - 

M87 (R<1';L01) 23.44±°°g - 23.52t°° g 1 - 23.30±°;°s 

NGC 4552 (l'<ii<10') 23.56 ± 0.20 1.34 ±0.12 23.75 ±0.35 1.40 ±0.18 23.33 ±0.16 1.09 ±0.10 

NGC 4552 (1'<K<5') 23.48 ± 0.22 1.42 ±0.14 23.85 ± 0.44 1.59 ±0.24 23.52 ± 0.28 1.31 ±0.18 

NGC 4552 (_R<1';KW01) 23.54 ±0.18 1.3 (fixed) - 

NGC 4552 (_R<1';L01) 23.19^'^ - 23.52^° ^ - 22.91±g ^| 



Table 2. Parameters of Gaussians fitted to GCLFs. In the outer region of NGC 4552, no Gaussian fits 
were attempted due to the poor statistics. The GCLF parameters in the core region (R < 1') of M87 and 
NGC 4552 taken from the literature are also shown for comparison: Kundu et al. (1999; K99), Kundu & 
Whitmore (2001; KW01), and Larsen et al. (2001; L01). L01 fitted t 5 functions to the GCLFs and hence we 
only indicate the Vro s - 




Figure 15. Histograms show GCLFs obtained only using GCs at distances $J 10' from host galaxy centre. Dotted line indicates 
incompleteness-corrected GCLF. Solid and dashed lines show GCLFs after subtracting control field populations based on the HDF-N 
and LH data, respectively. A Gaussian fitted to the GCLF is overplotted by a solid line. Shaded region indicates the magnitude range 
where the completeness is lower than 50 %. 



Vto = 23.56 ±0.20 and a = 1.34 ±0.12 mag. In Fig.EU the 
GCLFs of the red and blue GC subpopulations are shown 
with fitted Gaussians. Again the shaded region indicates the 
magnitude range where the completeness calculated by us- 
ing all the GCs (i.e., blue + red) is lower than 50 %. This 
limiting magnitude does not change a lot if only the red or 
blue GC subpopulation is considered since our B- and I-band 
images are deep enough not to miss a significant number of 
red or blue globular clusters detected on the V-band im- 
age. The fitted Gaussians suggest that the GCLF depends 
on subpopulation; the Vto of the red GC subpopulation is 
~ 0.5 mag and 0.4 mag fainter than that of the blue one 
for M87 and NGC 4552, respectively, although the results 
for the NGC 4552 GCs are less significant due to the large 
errors. Larsen et al. (2001) present similar results using the 
HST/WFPC2 data, although their Wo values tend to be 
brighter than those from other studies. This difference of 
Vto between the GC subpopulations is perhaps because of 
a metallicity difference (Ashman, Conti & Zepf 1995; El- 



son & Santiago 1996; Jordan et al. 2002). Vto and a of the 
GCLFs are summarized in Table [5] 

The GCs within 10' of M87 are further divided into two 
samples at a boundary of 5' and the GCLFs of all GCs, red 
GCs, and blue GCs in the inner and outer regions are pre- 
sented in Fig. 1171 Gaussians are also fitted to these GCLFs 
and the Vto and a are summarized in Table [2] This indi- 
cates that the shape of GCLF is not significantly different 
between the inner region (1' ^ R ^ 5') and the outer re- 
gion (5' ^ R ^ 10') for red GCs or blue GCs. Furthermore, 
the GCLF shape for all GCs is consistent with that in the 
core region 1') obtained by Kundu et al. (1999) using 
HST/WFPC2: Vto = 23.67 ± 0.07 mag (a = 1.39 ± 0.06). 
This suggests that the GCLF shape of the M87 GCs is not 
a strong function of distance from the host galaxy. We note 
that the fainter part of the GCLF for all GCs tends to be 
more deficient at the larger distance; this is due to the lower 
contribution of the red GC subpopulation in the outer re- 
gion. Any radial dependence of the GCLF is unclear for the 
GC population around NGC 4552 because the number of 
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Figure 16. The GCLF of all GCs, red GCs, and blue GCs is indicated by a black solid line, dotted line, and grey solid line, respectively. 
Incompleteness was corrected and contamination was subtracted based on the HDF-N data. Gaussians fitted to the GCLFs are also 
plotted. 




Figure 17. Same as Fig. 1161 but the GCLFs in the region of 1' ^ R ^ 5' around M87 are compared with those in the region of 
5' ^ R ^ 10'. 



GCs is not substantial enough especially in the outer region 
(5' s; 7? ^ 10'). Nevertheless, the Vto of the GCLF obtained 
at l' ^ R ^ 10' is consistent with that in the core region ob- 
tained by Kundu & Whitmore (2001) using HST/WFPC2 
who found Vto = 23.47 ± 0.14 mag with a fixed value of 
a — 1.3, again suggesting that the GCLF shape does not 
depend significantly on distance from the host galaxy. 

4.2 GC Specific Frequency 

Using the Vto and a of the Gaussians fitted to the GCLFs 
presented in Fig. 1151 we can estimate the total number of 
GCs associated with M87 or NGC 4552 and calculate a GC 
specific frequency (Sn)- We first integrate a GCLF obtained 
in an armulus (1' width) centered on the host galaxy down 
to the magnitude where the completeness becomes 50 %. 
This number is then multiplied by a correction factor to in- 
clude the fainter GCs based on the Gaussian GCLF. This 
calculation is repeated out to a certain distance from the 
galaxy centre. We note that it is ideal to fit a Gaussian 



to a GCLF obtained in each annulus but this is difficult in 
practice, especially at large distances because the number of 
GCs is not substantial and statistical errors are significant. 
We therefore use the fitted Gaussian in the inner region 
10') of M87 or NGC 4552 independently of the distance 
from the host galaxy. In order to estimate the number of 
GCs at distances smaller than ^ 1', where the brightness 
of host galaxy halo light exceeds the linearity regime of the 
CCD, the radial profile of GC surface density outside the 
core region is fitted with a de Vaucouleurs law profile (see 
Paper II for details) and the fitted formula is extrapolated 
towards the galaxy centre. The estimated total number of 
M87 GCs is 12000 ± 800 within 25' from the galaxy centre 
(cf. 13200 ± 1500 by HHM). For NGC 4552, the total num- 
ber of GCs is estimated to be 1400 ± 170 within 10'. Note 
that the fraction of GCs at R ^ 1' estimated from the ex- 
trapolation is 8.5 % for M87 GCs and 17 % for NGC 4552 
GCs. If we adopt Mv = —22.46 mag as the V-band absolute 
magnitude of M87, which is also adopted by HHM (a differ- 
ence in adopted distance modulus of 0.03 mag is corrected), 
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the Sn value of M87 is calculated to be 12.5 ± 0.8, which 
is only slightly smaller than 14.1 ± 1.6 obtained by HHM. 
The T/-band luminosity of NGC 4552 is obtained from our 
data by fitting a de Vaucouleurs law to the V-band surface 
brightness profile and integrating it out to the same distance 
(10'); M v = -21.12 mag. The S N value of NGC 4552 is then 
estimated to be 5.0 ± 0.6. 4 



5 SUMMARY 

We have performed a wide-field imaging survey of the glob- 
ular cluster (GC) populations around M87 with Suprime- 
Cam on the 8.2m Subaru telescope. A 2° x 0°5 (560 kpc x 
140 kpc) field extending from M87 to the east was observed 
through the BVI filters. In addition to this unprecedented 
large survey area, our data analysis has been optimized to 
study the statistical properties of GCs as follows: 

• GC candidates are isolated not only with an extended 
source cut but also with a colour cut, where only unresolved 
objects falling on a specific region of the B — V and V — I 
colour-colour diagram are accepted. The colour criterion is 
defined so as to include almost all the Galactic GCs and 
to avoid foreground stars and background galaxies. This is 
expected to efficiently isolate bona-fide GCs from other un- 
resolved objects on our imaging data. 

• In order to assess foreground and background contami- 
nation which needs to be statistically subtracted, we analyze 
the imaging data on the HDF-N field and the Lockman Hole 
field as control field data. These fields cover reasonably wide 
sky areas (~ 30' x 30') and are compatible with the data of 
the M87 fields in terms of the filter set (BVI), limiting mag- 
nitudes, and image qualities. We therefore extract contami- 
nating populations using identical criteria to those adopted 
in the M87 fields, minimizing the possibility of introducing 
any systematic errors into the subtractive correction. 

In this paper, we have investigated the luminosity func- 
tion and global specific frequency (Sn) of GC candidates 
surrounding M87 or NGC 4552. The V'-band GC luminos- 
ity functions (GCLFs) were obtained in the inner regions 
of M87 and NGC 4552 at distances < 10' from the galaxy 
centres. By fitting Gaussians to the GCLFs, the turnover 
magnitude is estimated to be 23.62 ±0.06 mag for M87 GCs 
and 23.56 ± 0.20 mag for NGC 4552 GCs. The GCLF ap- 
pears to depend on GC colour; the turnover magnitude in 
the GCLF of the red GC subpopulation (V — I > 1.1) is 
~ 0.5 mag and 0.4 mag fainter than that of the blue GC 
subpopulation (V - I < 1.1) for the M87 GCs and NGC 
4552 GCs, respectively. 

For the M87 GCs, the GCLFs at 1' < R < 5' were 
compared with those at 5' ^ R 5j 10' but no obvious trend 
with radius was found in the shape of the GCLF for either 
the red or blue subpopulations. The global Sn of M87 GCs 

4 We calculated the total number of GCs and Sn within 10' 
for NGC 4552 to avoid possible contributions of M87 GCs and 
intergalactic GCs outside of this radius (see Paper II). If we fit a 
de Vaucouleurs law to the GC surface density profile within 10' 
from the NGC 4552 centre and integrate it out to 25' as done 
for M87 GCs, the number of GCs and Sn are estimated to be 
2000 ± 660 and 6.2 ± 2.1, respectively. 



and NGC 4552 GCs is estimated to be 12.5 ± 0.8 within 25' 
and 5.0 ± 0.6 within 10', respectively. 
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